Non-Photonic Electron pt Distributions and Correlations of Electrons from B and D 

Meson Decays with Charged Hadrons 
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We compare the non-photonic electron and D° meson pr distributions measured at RHIC with the 
PYTHIA Monte Carlo event generator in p+p collisions at y^ijvjv = 200 GeV. A delta fragmentation 
function much harder than the Peterson function, consistent with recombination scheme for charm 
meson formation, is needed to match the experimental data. Attempts to fit the non-photonic 
electrons at high pT show large uncertainties and the B meson semi-leptonic decays may not be 
dominant for electron pt up to 8 GeV/c. Correlations of non-photonic electrons with charged 
hadrons are studied. We propose an experimental method to quantitatively determine the relative 
contributions of D and B meson semi-leptonic decays to the non-photonic electrons. 
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Heavy quarks are believed to be produced through 
initial parton-parton, mostly gluon-gluon, scatterings in 
nuclear collisions at the Relativistic Heavy Ion Collider 
(RHIC) energies. Theoretical calculations of heavy quark 
production within the perturbative Quantum Chromo- 
Dynamics (pQCD) framework are considered more reli- 
able because the heavy quark mass sets a natural scale for 
the pQCD. Charm quark production is sensitive to the 
incoming parton flux for the initial conditions of nuclear 
collisions The transport dynamics of the heavy 

quarks in nuclear medium such as flow and energy 
loss 00 can probe QCD properties of the dense matter 
created in nucleus-nucleus collisions. Therefore, heavy 
quark measurements provide unique insights into QCD 
properties of the new state of matter produced in nucleus- 
nucleus collisions. 

The heavy quark production in p-l-p or p-l-A collisions 
provides a reference for heavy meson formation and for 
nuclear modification factors of heavy quarks in the nu- 
clear medium. The Peterson function (61 
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has often been used to describe the charm fragmenta- 
tion function, where the parameter e « 0.05 (default 
in PYTHIA 0) is in reasonable agreement with the 
results from fits to charm production in e+e^ and 773 
collisions 0. However, in charm hadroproduction, it 
was observed that the c-quark pT distributions of next- 
to-leading-order (NLO) pQCD calculations agree well 
with the measured open charm px spectrum Q , indicat- 
ing that a much harder fragmentation function peaked 
at z « 1 in Eq. ^ is needed in charm hadroproduc- 
tion. A more detailed discussion of this observation 
can be found in [lol|[ll| . Coalesence or recom- 

bination (U] |0| 01 models have also been proposed for 
charmed meson formation by combining a charm quark 
with a light up or down quark, presumably of soft px |l7l | . 



Thus the charmed hadron pT would coincide with the 
bare charm quark pT distribution in this hadronization 
scheme. These various hadron formation schemes can 
lead to significantly different charmed meson pT distribu- 
tion when interpreting non-photonic electron px spectra 
from experimental measurements. 

Recently STAR and PHENIX collaborations have ob- 
served a suppression of non-photonic electrons much 
larger than predicted in the pt ~ 4 — 8 GeV/c re- 
gion 10 01 . This observation challenges theoretical pre- 
dictions based on energy loss via induced gluon radia- 
tion for heavy quarks . The fact that the relative 
contributions of D and B meson semi-leptonic decays to 
the experimental non-photonic electron spectrum is cur- 
rently not measured at RHIC severely restricts the under- 
standing of the heavy quark energy loss. A nonzero non- 
photonic electron V2 has been measured in the pr < 2.0 
GeV/c region, while at higher the V2 is observed to de- 
crease with Pt [12 ■ Quantitative understanding of this 
feature in heavy quark measurements requires relative 
charm and bottom contributions to non-photonic elec- 
trons. 

In this paper we evaluate the px distribution of 
D mesons from PYTHIA v6.22 and compare the 
PYTHIA results with the STAR measurements |2^[2^. 
The charm quark fragmentation function will be mod- 
ified from the default Peterson function and the other 
PYTHIA parameters are tuned in order to describe the 
experimental D meson data as well as the STAR non- 
photonic electron spectra [3] ■ We study the correlations 
between non-photonic electrons and charged hadrons in 
order to guide the experimental investigations of the rela- 
tive contributions to the non-photonic electron spectrum 
at high Pt from B and D meson decays. 

Fig. ^ shows charm quark and meson spectra 
from PYTHIA calculations and NLO pQCD predictions 
for charm quark spectra [i^ together with the STAR 
D'^ spectrum [i^, where the measured data points 
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FIG. 1: (color online) Charm quark (left) and (right) 
spectra from PYTHIA calculations compared with next-to- 
leading-order pQCD predictions for charm quark spectra and 
the STAR measured spectrum from d-|-Au collisions scaled 
by Ntin = 7.5. 



from d-l-Au collisions at y^SNN = 200 GeV have been 
scaled by iVbi„ = 7.5 - the number of binary collisions. 
The PYTHIA spectra have been scaled to the measured 
dN/dy of 0.028 ±0.004(siai.)±0.008(sysi.) for D" at 
mid-rapidity, where the scaling factors are 1.75 for 1)° 
meson spectra and 1.25 for charm quark spectra, respec- 
tively. The theoretical curve with the MRST HO PDF 
has been normalized to the measured total cc cross sec- 
tion (1.3 ± 0.2 ± OAmb) ^ by a factor of 3.4. 

The stars in the left panel of Fig. ^ depict the 
charm quark px distribution from a PYTHIA calcula- 
tion with the following parameters (set I) from refer- 
ence 10: PARP(67) = 1 (factor multiplied to Q^) ^ 
< kt >= 1.5GeV/c, TOc = 1.25Gey/c2, Kfactor = 3.5, 
MSTP(33) = 1 (inclusion of K factors), MSTP(32) ^ 
4 (Q2 scale) and CTEQ5L PDF. We further tuned the 
value of PARP(67) to 4, which enhances c-quark produc- 
tion probability through gluon splitting and is introduced 
to take into account higher order effects in the pQCD cal- 
culation 27]. The results are shown as triangles in the 
left panel of Fig. ^ We found that this change mainly 
affects c-quark production at high px and can effectively 
reproduce the NLO pQCD calculation. We will refer to 
this set of parameters as parameter set 11 in the rest of 
the paper. 

We compared the charm quark spectra from PYTHIA 
calculations to the STAR published spectrum 1231 
together with the STAR preliminary D* spectrum |24| . 
The combination of STAR measured D* and D° me- 
son spectra into a single D meson spectrum is ham- 
pered by the large uncertainties in the ratio D* /D^ — 
0.4 ± 0.09 ±0.13 which is experimentally not well 



FIG. 2: (color online) Electron spectra from PYTHIA calcu- 
lations with the modified Peterson function for charm quarks 
and beauty quarks compared with background subtracted sin- 
gle electron spectrum measured by STAR from p-f p collisions 
and the prediction for theoretical uncertainty electron band. 
Left panel: The spectra are from parameter set II. Right 
panel: The spectra are from parameter set I. 

known at RHIC. The STAR D meson spectrum covers 
a limit pT range and is not well constrained at high px- 
Without further tuning other PYTHIA parameters, we 
find that the generated bare charm quark spectra from 
PYTHIA calculations using parameter set I and param- 
eter set 11 approximately match the STAR D meson pT 
distribution. As demonstrated in the left panel of Fig.^ 
the PYTHIA calculation with parameter set 11 also yields 
a charm quark distribution similar to the NLO pQCD 
calculations (251] (22|], which coincides with the STAR D 
meson pT spectrum as well. It is not the purpose of 
this paper to show that the PYTHIA calculations with 
these two parameter sets and NLO pQCD calculations 
are equivalent in physics contents. 

In the right panel of Fig. Q] the D° pt distribution 
from PYTHIA calculations using the default Peterson 
fragmentation function is shown as open circles for pa- 
rameter set 11 and crosses for parameter set 1. The default 
Peterson function refers to the value of the parameter e in 
Peterson function being 0.05 for charm quarks and 0.005 
for beauty quarks. The default Peterson fragmentation 
for charm quarks is too soft to reproduce the measured 
spectrum together with D* spectrum ^ij- We 
modified the value of the parameter e to 10~^ for both 
charm and beauty quarks. In this case the fragmentation 
function is nearly i5(l — 2). The results are shown as stars 
for parameter set 1 and triangles for parameter set 11 in 
the right panel of Fig. ^ The PYTHIA calculations using 
the modified Peterson fragmentation function {e = 10"^) 
with parameter set 1 and set 11 can reasonably reproduce 
the measured D meson px distribution. 
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FIG. 3: (color online) Atp distributions between non-photonic 
electrons and charged hadrons with diflterent electron trigger 
Pt cuts and associated hadron pr > 0.1 GeV/c. Solid lines 
show electrons from B meson decays. Dashed lines show elec- 
trons from D meson decays. 



FIG. 4: (color online) Summed pT distributions of charged 
hadrons around triggered non-photonic electrons with differ- 
ent trigger px cuts. Solid lines show electrons from B meson 
decays. Dashed lines show electrons from D meson decays. 



While the kx broadening can make the charm pT dis- 
tribution harder at low beam energies the intrinsic 
kr has little effect on the pr distribution at RHIC ener- 
gies ^25]. A harder fragmentation function is needed for 
the hadronization of charm quarks if the pQCD calcula- 
tion is to describe the measured STAR data. 

The STAR independent measurements of the recon- 
structed D° meson and single electrons from heavy quark 
semi-leptonic decays measured with TOF and TPC are 
consistent The electron measurement there only 

covers up to < 4 GeV/c and has no sensitivity to the 
B meson contribution. We also checked the consistency 
between D meson data and non-photonic single electron 
data in p-l-p collisions within our PYTHIA calculations. 
Fig. El shows the electron spectra from PYTHIA calcula- 
tions using the modified Peterson fragmentation function 
for charm quarks and beauty quarks with parameter set 
II in the left panel and parameter set I in the right panel. 
The parameters for beauty quarks in set II and set I are 
all the same as for charm quarks except rUb = i.SGeV/c^. 
The PYTHIA spectra of electrons from charm meson 
decays are scaled by the same factor used to scale the 
PYTHIA D° spectra to the measured dN/dy{D°) at mid- 
rapidity. The electron spectra from beauty meson decays 
are normalized by the ratio of abi/ucc based on the NLO 
pQCD calculation [i^. The band corresponds to the 
theoretical uncertainty of this ratio (0.45% - 0.60%) [j^ . 
where the theoretical error may be an underestimate ac- 
cording to recent calculations We used the value 
at the center of this range to calculate the sum of the 
electrons. Fig.|21also shows the comparison to the STAR 
measured non-photonic single electron data in p-f p colli- 
sions fl^ and the FONLL calculation for the theoretical 
uncertainty band of the electron spectrum from charm 



and bottom in p-l-p collisions |23|. The FONLL predic- 
tion of electron spectrum gives a fair description of the 
shape of the measured spectra. But there is a discrepancy 
in the overall scale and the FNOLL calculation is signifi- 
cant below the STAR data at high px- More discussions 
can be found in 18]. The PYTHIA calculations with pa- 
rameter set II and modified heavy quark fragmentation 
function can simultaneously describe the STAR direct 
D meson pr distribution l2il2l and the STAR EMC 
non-photonic electron data 18]. This modified scheme 
indicates that the contribution of electrons from B me- 
son decays is not dominant for the measured px region 
up to 8 GeV/c assuming (7^1! (Jcc ~ 0.45% — 0.60% based 
on NLO pQCD calculation. The measurement of elec- 
tron Pt distribution alone has a reduced sensitivity to the 
Pt distribution of D mesons as shown in reference '30|. 
Within the statistical and systematic errors of the elec- 
tron data, the PYTHIA calculation with parameter set 
I and the modified Peterson function for both c and h 
quarks can yield an electron pT distribution similar to 
the measurement |0 . In this case, as shown in the bot- 
tom panel of Fig.|21 electrons from h quark decays are not 
the dominant source for the pT region up to 6 GeV/c. It 
is critical that the B and D meson decay contributions to 
non-photonic electrons to be separated experimentally. 

We have studied the azimuthal correlations between 
heavy quark semi-leptonic decay electrons and inclusive 
charged hadrons. Fig. shows the At/? distributions 
between non-photonic electrons and inclusive charged 
hadrons, with different electron trigger pT cuts and the 
associated hadron > 0.1 GeV/c. The distributions are 
scaled by the number of electron triggers. These plots are 
from PYTHIA calculation with parameter set II and the 
modified heavy quark fragmentation function. The solid 
lines in Fig.|31are for electrons from B meson decays and 
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TABLE L Fractions of B meson decay contributions from 
electron spectra and Aip distribution fitting in PYTHIA cal- 
culations with parameter set 11 and modified Peterson frag- 
mentation function for charm and beauty quarks. 



Pritrig) (GeV/c) 


From Spectra (%) 


From Fitting (%) 


2.5 - 3.5 


12.64 ±0.04 


12.64 ± 1.37 


3.5 - 4.5 


21.65 ±0.12 


21.65 ±2.55 


4.5 - 5.5 


30.50 ± 0.26 


30.50 ±4.76 


5.5 - 6.5 


36.66 ± 0.48 


36.66 ± 7.87 


6.5 - 8.5 


42.24 ±0.71 


42.24 ± 11.10 


8.5 - 10.5 


49.82 ± 1.72 


49.82 ±25.11 



the dashed lines are for electrons from D decays. As 
demonstrated in Fig. |31 there is a significant difference 
between B and D meson decays in the near-side corre- 
lations. The width of near-side peak for electrons from 
D decays is much narrower than those for the B decays. 
The wide width from B meson decays is due to the larger 
energy release {Q value) in the B meson semi-leptonic 
decays leading to a broad angular correlation between 
daughter hadrons and electrons. An electron at high px, 
if it is from B meson decays, the B meson does not have 
to be at high transverse momentum because the electron 
can get large momentum from the 6-quark mass. In the 
case of D meson decays, the D meson needs to have a 
large momentum in order to boost the daughter electron 
to a high pt- In experimental effort we cannot differen- 
tiate heavy quark decay hadrons from others, we used 
inclusive hadrons in our correlation study. We found the 
difference in the near-side correlations between D decays 
and B decays is largely due to the decay kinematics, not 
the production dynamics. This difference can help us es- 
timate the relative B and D contributions to the yields 
for non-photonic electrons. 

The bottom contribution to the non-photonic electrons 
can be determined directly from the electron spectra in 
PYTHIA. The results are given in the second column of 
TableQ] However, this way to determine the bottom con- 
tribution is not experimentally feasible. In order to ex- 
perimentally determine the bottom contribution fraction, 
we use the Aip distributions for B decays and D decays 
to fit the A(p distribution for PYTHIA inclusive case (ei- 
ther B decays or D decays), and let the B contribution 
fraction as a parameter. The fraction is determined by 
the minimum value of x^- And the fitting error is deter- 
mined by one a shift of /ndf from the minimum value. 
The fitting error will reduce by increasing the statistics. 
The results are shown in the third column of Table 
The results are consistent with those directly from the 
electron spectra. This indicates that the method we pro- 
posed is self-consistent. This approach can be used to 
experimentally estimate the B and D contributions to 
the non-photonic electrons. 



TABLE II: Fractions of B meson decay contributions from 
electron spectra and summed pr histogram fitting in PYTHIA 
calculations with parameter set II and modified Peterson frag- 
mentation function for charm and beauty quarks. 



PT{trig) (GeV/c) 


From Spectra (%) 


From Fitting (%) 


2.5 - 3.5 


8.25 ±0.05 


8.25 ±1.50 


3.5 - 4.5 


14.94 ±0.13 


14.94 ± 2.54 


4.5 - 5.5 


22.60 ±0.29 


22.60 ±4.13 



We further studied the particle production within a 
cone around triggered high px electrons from heavy quark 
decays. We focused on the scalar summed pr distri- 
butions of inclusive charged hadrons in the cone {pt 
refers to the transverse momentum in the laboratory 
frame). Here the cone is defined by jry^ — fle\ < 0.35 
and \iph — f>e\ < 0.35 (77 is pseudorapidity and ip is az- 
imuthal angle). The summed px distributions of inclu- 
sive charged hadrons in three triggered electron px ranges 
are shown in Fig. 01 The distributions are scaled to unit. 
The dashed lines are for D decays and the solid lines are 
for B decays. These plots are from PYTHIA calcula- 
tion with parameter set II and the modified heavy quark 
fragmentation function. We also can see that there is 
a significant difference between B decays and D decays. 
The summed pt distributions for D meson decays are 
much wider than those for B meson decays. This dif- 
ference can also be used to distinguish B and D decay 
contributions. We use the summed pr histograms from 
B decays and D decays to fit the summed pr histogram 
from PYTHIA inclusive case to determine the B decay 
contribution. The results are shown in Table ITTl The re- 
sults are consistent with those directly from the electron 
spectra. The ratios for the same trigger pT ranges are 
different between Table IhI and Table It is because we 
removed those electrons which have no hadrons in the 
cone around them when we calculated the bottom con- 
tribution from the electron spectra for Table ^ This 
removes different fractions from B and D decays, which 
has to be corrected by simulations. This was done to 
make the results directly from the electron spectra com- 
parable to the results from hadron summed pr histogram 
fitting. It will be more valuable for the small acceptance 
experiments to investigate the B decay contribution us- 
ing summed px histogram fitting. 

In conclusion, we found that in order to match the 
STAR measured pr shape of D mesons from d+Au col- 
lisions at RHIC using the PYTHIA Monte Carlo gener- 
ator, a harder charm quark fragmentation function must 
be used. Since the D meson production at high pT is 
not well constrained by the STAR data and the mea- 
surement of electrons alone has a reduced sensitivity on 
the Pt shape of D mesons, PYTHIA calculations with 
parameter set I and set II can simultaneously describe 
the Pt distributions of D mesons and non-photonic elec- 
trons from semi-leptonic decays of heavy quarks. These 
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calculations predict that the beauty quark decay elec- 
trons are not a dominant contribution over the entire pt 
region up to 6-8 GeV/c depending on parameters. The 
relative contributions to non-photonic electrons from B 
and D decays have to be determined experimentally. We 
studied the correlations of non-photonic electrons and in- 
clusive charged hadrons, which can distinguish between 
D and B decay contributions to the non-photonic elec- 
trons due to their different decay kinematics. We have 
proposed an experimental method to estimate the rela- 
tive D and B decay contributions. 
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